Quantum interference and evolution of entanglement in a system of three-level atoms 
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We consider a pair of three-level atoms interacting with the vacuum. The process of disentanglement due 
to spontaneous emission and the role of quantum interference between principal transitions in this process, are 
analysed. We show that the presence of interference can slow down disentanglement. In the limit of maximal 
interference, some part of initial entanglement can survive. 
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I. INTRODUCTION 

The process of decoherence and degradation of entangle- 
ment in open quantum systems was studied mainly in the case 
of pair of two-level atoms interacting with vacuum or thermal 
noises. This system is simple to study, since the amount of 
its entanglement can be quantified by concurrence 1 1 ] . This 
quantity discriminates between separable and entangled states 
and is analytically computable. The dynamics of the system 
given in the Markovian regime by the master equation, can be 
used to obtain time evolution of concurrence and to analyse 
the influence of the process of decoherence on entanglement. 
Let us mention two particular results of such studies. The pro- 
cess of disentanglement was considered in the case of a pair 
of atoms, each interacting with its own reservoir at zero tem- 
perature. It turns out that concurrence of initially entangled 
state can vanish asymptotically or in finite time, depending 
on initial state |2, 3, 4]. This behaviour should be compared 
with asymptotic decoherence of any initial state. On the other 
hand, individual atoms located inside two independent envi- 
ronments at infinite temperatures, always disentangle at finite 
times 01 . 

In the present paper, we extend this kind of studies to the 
case of pair of three-level atoms. The analysis is much more 
involved since there is no simple necessary and sufficient cri- 
terion of entanglement for a pair of c/-level systems with t/ > 3. 
Peres' separability criterion |6] only shows that states which 
are not positive after partial transposition (NPPT states) are 
entangled. But there can exist entangled states which are pos- 
itive after this operation |7] (PPT states). Using this criterion 
we at most can study the evolution of NPPT states and ask 
when they become PPT states, but in the models considered 
in the present paper we can directly observe when stationary 
asymptotic states are separable. 

Interesting example of dissipative dynamics comes from 
the analysis of three-level V-type atomic system where spon- 
taneous emission may take place from two excited levels to 
the ground state and direct transition between excited levels is 
not allowed. However the indirect coupling between excited 
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States can appear due to interaction with the vacuum (quan- 
tum interference) |8]. This interference results from the fol- 
lowing mechanism: spontaneous emission from one transition 
modifies spontaneous emission from the other transition j^l- 
There were many studies on the effect of quantum interfer- 
ence on various processes including: resonance fluorescence 
fioll . quantum jumps the presence of ultranarrow spec- 
tral lines fl2ll or amplification without population inversion 

In our research we consider a pair of such three-level 
systems and study another aspect of quantum interference, 
namely its influence on the evolution of entangled states. In 
the case of distant atoms we expect that dissipation causes 
disentanglement, but the rate of this process may depend on 
interference. It is indeed true, as shows analysis of the models 
considered in the paper. We demonstrate it for some classes 
of pure and mixed states by proving that the measure of entan- 
glement vanishes asymptotically, but the larger is the effect of 
interference, the slower is the process of disentanglement. In 
the limit of maximal interference we obtain very interesting 
phenomenon: some part of initial entanglement can survive 
and we obtain asymptotic states with non-zero entanglement 
despite of of the inevitable process of decoherence. 

In this paper we study two models of three-level systems 
with quantum interference between transitions. The simpler, 
proposed in Ref. 1 9], applied to the pair of atoms can be solved 
analytically. This model simulates many aspects of V-type 
three-level atom where interference effect was discovered 1 8|] . 
We study the dynamics of a pair of such atoms numerically 
and show that both models predict similar behaviour of entan- 
glement. 



II. THREE-LEVEL ATOM AND QUANTUM 
INTERFERENCE 

Consider a three-level atom in the V configuration. The 
atom has two nondegenerate excited states |1), |2) and the 
ground state |3). Assume that excited levels can decay to the 
ground state by spontaneous emission, whereas a direct tran- 
sition between excited levels is not allowed. If the dipole mo- 
ments of these two transitions are parallel, then indirect cou- 
pling between states |1) and |2) can appear due to interaction 
with the vacuum (quantum interference between transitions 
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|1)^|3) and |2) ^ |3)) (FIG. 1). 
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FIG. 2: Three-level atom in the V configuration (system II) 



FIG. 1: Three-level atom in the V configuration (system I) 

Time evolution of such system (which we call system I) is 
given by the master equation 1 9J 



^ = -/[//, p]+L/P 



(II. 1) 



where the damping term is 
1 



^/P= 2"^^ (203lP0l3-0iip-p0li) 

+ ^ 72 (2032PO23 - a22P - pa22) 

1 , 

+ 2 yi2(203lP(J23 - O21P - P02l) 



(11.2) 
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yi2(2032POi3 - O12P - PO12) 



In this equation, Ojk is the transition operator from \k) to | /) 
and Yij Y2 are spontaneous emission constants of |1) and |2) 
to the ground level |3). In addition 



Y12 



(11.3) 



gives cross damping term between |1) |3) and |2) |3). 
The parameter p/ represents the mutual orientation of transi- 
tion dipole moments: when p/ = 1, quantum interference is 
maximal and for P/ = it vanishes. Since we are mainly in- 
terested in the evolution of initial states due to the spontaneous 
emission, the form of the Hamiltonian H is not discussed (see 
1^ for details). 

It is known that in atoms used in atomic spectroscopy, the 
transition dipole moments are usually perpendicular, so the 
system (I) described by equation ( III.2> is difficult to realize. 
However, as was shown in Ref. |9], the effects of quantum 
interference can be duplicated to a large degree by the three- 
level system (II) (FIG. 2) with exited states |e), \u) and the 
ground state \g), evolving according to the master equation 
with damping operator 



+ l-ju (20gi,po„g - OuuP - pa„„) 



(11.4) 



In the system (II), the dipole moments can be perpendicular 
and the measure of quantum interference is given by 



P// 



(11.5) 



So we can expect maximal effects of quantum interference in 
the system (II) when the level \u) is metastable. 



III. ENTANGLED PAIR OF THREE-LEVEL ATOMS 

Now we consider a pair of three-level systems in the V 
configuration. In the context of evolution given by master 
equation generated by (III.2> or MIA\ . we address the follow- 
ing question: how the dynamics of compound system of two 
atoms prepared initially in the entangled state is influenced by 
the presence of quantum interference? In particular, we may 
look at the time evolution of the appropriate measure of en- 
tanglement. Since damping causes that pure states evolve into 
mixed states, we need effective measure of mixed - state en- 
tanglement. As such measure one usually takes entanglement 
of formation £f(p) O, but in practice it is not known how 
to compute this measure for mixed states of pairs of d - level 
systems in the case when d > 2. A computable measure of 
entanglement was proposed in Ref. lITsll . It is based on the 
trace norm of the partial transposition p'^ of the state p. From 
the Peres' criterion of separability |6], it follows that if p^'* is 
not positive, then p is not separable. So one defines negativity 
of the state p as 



A^(p) = 



lp^-lli-1 



(III.l) 



A^(p) equals to the absolute value of the sum of negative eigen- 
values of p^'* and is an entanglement monotone L151 . but it can 
not detect bound entangled states 0]. 

The compound system of two atoms is defined on the space 
§5 and the matrix elements of the state p will be consid- 
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ered with respect to the basis 

{ei (g) ei, ei (g) 62, ei (g) 63, 62 g) ei, 62 (g 62, 62 63, 
63 (g) ei, 63 (g 62, 63 g) 63} 

(III.2) 

where {ei, 62, 63} is the canonical basis of C^. In our dis- 
cussion of entanglement evolution we focus on two classes of 
initial states. The first one consist of pure states of the form 

^9 (p = cos9sin(pei g) ei + sin0sin(pe2 g) 62 + cos 963 g) 63 

(III.3) 

where 9, 9 G [0,7t/2]. The states (IIII.3> have negativity 

N{^Q ^ ) = cos 9 sin (p cos cp + cos 9 sin sin'^ (p + sin 9 cos cp sin (p 

(111.4) 

In the special case 9 = 7t/4, (p = arccos we obtain the state 



where 



V3 



(ei(gei 



- 62 g) 62 + 63 (g 63 j 



(III.5) 



with maximal negativity A^CPmax) = 1- The second class of so 
called isotropic mixed states Iil6l . is defined as follows: let 

PO = y 

where Ig is the identity matrix in C^, be the maximally mixed 
state of two three - level systems, then 



W = (l-p)pO + p|^max)(^n 



pe[o,i] (111.6) 



interpolate between maximally entangled and maximally 
mixed states. Notice that states (IIII.6> are the natural gen- 
eralizations of Werner states 1 17] of two-qubit systems. One 
can check that negativity of ( IIII.6t equals 



1. 



Moreover, we identify the states |e), \u) and \g) with vectors 
61, 62 and 63 respectively. The master equation ( IIV.H can 
be used to obtain the equations for matrix elements of any 
density matrix of compound system, with respect to the basis 
( IIII.2> . The calculations are tedious but elementary, so we are 
not discussing details. As a result, we obtain the following 
properties of dynamics given by dlV.H : 

1. for jii > 0, the dynamics brings all initial states into 
unique asymptotic state \g) g) \g), 

2. in the case of maximal interference (y„ = 0), the dynam- 
ics is not ergodic, and asymptotic stationary states Pas 
depend on initial conditions. 

The case of maximal interference will be studied in the next 
section, now we discuss the evolution of entanglement in the 
first case. Since the asymptotic state is separable, negativity 
of any initially entangled state goes to zero with time. To see 
some details, consider initial states of the form 



P51 



VP91 



Pl5 
P55 
P95 



Pl9\ 

P59 

P99y 



(IV.3) 



Notice that pure states JIII.3> are particular examples of ( IIV.3> . 
The solution p(f) of master equation dlV. 1> for initial condi- 
tion (IIV.3> has matrix elements 



N{W) : 



0, 



■(4p-l), p>i 



(111.7) 



IV. TIME EVOLUTION OF NEGATIVITY 

Consider first the system of two distant three - level atoms 
A and B , both of type (II). This case is simpler to analyse and 
can be solved exactly. Since the atoms are independent and 
do not interact, the dissipative part of dynamics is given by 
the master equation 



dp _ jAB^ 



(IV. 1) 



with 



4^P-^Y.(2<P< 


-a^,p- 


~pO 


+ \lu (2<,p< 


-<^,uP- 


-P<<u) 


+ \le{2(5%P<^%' 


'^eeP~ 


-pofj 


+ \lu (2of„paf^, 


-oIp- 


-P<u) 



(IV.2) 



pii(0 


= e 


'^^'Pil 


Pl5(f) 


— 




Pl9(f) 
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^''Pl9 


P33(0 


^{e 




P55(0 




'^^'P55 


P59(f) 
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^"'P59 


P66(r) 


^{e 


-V_e-2V)p55 


P77(r) 


^{e 


-V-e-2V)p„ 


P88(f) 


= {e 


-W_,-2V)p33 


P99(0 


= (1 






+ (1 





(IV.4) 



where for remaining matrix elements one can use hermiticity 
of p(f). Even in that case, the analytic formula for negativity 
as a function of time is rather involved, so we consider ex- 
plicit example. Take as initial state ^max which is maximally 
entangled, then 



(IV.5) 
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FIG. 3: Time-dependence of negativity for initial state "Pmax for dif- 
ferent values of a = Y«/Ye 
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FIG. 4: Time-dependence of negativity for initial state <III.3t with 
6 = Jr/8, (p = k/6 and different values of a = Y«/Ye 



Observe that (IIV.5t vanishes asymptotically at the rate de- 
pending on degree of interference: the larger is the effect 
of interference, the slower is the process of disentanglement 
(FIG. 3). Notice also that in the limit of maximal interfer- 
ence (y„ — > 0), some part of initial entanglement can survive. 
One can also check that other entangled states from the class 
jIII.3> behave similarly (see FIG. 4). When the initial state is 
mixed isotropic state W given by ( IIII.6> . then has matrix 
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FIG. 5: Time-dependence of negativity for initial state W with p - 
3/4 and different values of a = Ju/je 



elements 



W„(f)-^(l+2p)e-2V 
Wi9(f) = 



W22(f) 
W33(f) 
W44(f 

W59(f 



) = i(l_;,)e-(Tm«)' 



-Jet 



{l+2p)e-^'f"' 



W66{t) = ^{p-l)e-^''^+''"'>'- 

W7i{t)^W33{t) 
^88(0=^^66(0 



^{l+2p)e 



W99(f 



(IV.6) 



Explicit formula for negativity of W{t) is lengthy and in- 
volved, so we not reproduce it here. But it follows from the 
calculations, that the time dependence of negativity is similar 
to the case of pure states (FIG. 5). 



V. MAXIMAL INTERFERENCE AND ASYMPTOTIC 
ENTANGLEMENT 

The case of maximal interference in system (II) can be dis- 
cussed in details. Take any initial state p with matrix elements 
p,y. The dynamics given by the master equation (IIV.U has a 
remarkable property: in the case of maximal interference be- 
tween transitions \e) \g) and \u) \g), there are nontrivial 
stationary asymptotic states pas. By a direct calculations, one 
can show that Pas have the following non-vanishing matrix el- 
ements 



(Pas)55=P55, (Pas)56 = P56, 
(Pas)59 = P59, (Pas)68 = P68, 
(Pas)69 = P47 + P69, (Pas)88 = 
(Pas)89 = P23 + P89, (Pas)99 = 



(pas)58 = P58, 
(pas)66 = p44 + p66 

- P22 + P88 

- Pll +P33 + P77 + P99 

(V.l) 



It is very interesting that although two atoms are independent 
and do not interact, some of these asymptotic states can be 
entangled. Quantum interference causes that part of initial 
entanglement can survive the process of decoherence. For ex- 
ample, if initial states are pure states (IIII.3> . then Pas has only 
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four nonvanishing matrix elements 

(Pas)55 = sin^cpsin^e, (Pas)59 = COS cpsin (psinG, 
(Pas)99 = cos^(p + sin"(pcos^9 

and 

A^(Pas) — sin (p COS (p sin 9 
Comparing (|n^3 with ( IV. 3> we see that 

A^(pas) < A?(^e,<p) 

Consider now some special cases: 
1. Let 9 = 7t/4 then 



(V.2) 



(V.3) 



sm(p 

71 



(ei ei + 62® 62) + cos (pes (g) 63 



where < (p < k/2. The states have negativity 
N{^'(f) = isin(p(2\/2cos(p + sin(p) 



and 



A^(pa 



smcp cos(p 
V2 



which is always greater then zero. 

2. Let 9 = ±71/2, (p e (0,7t/2), then the states 

= ±sin(pe2§5e2 + cos (p 63 (g) 63 
are stable during the evolution, and 

A^(Pas) ^ N(^^) = coscp sincp 

3. Let 9 = 0, then the states 

"PjJ = sincpei +cos(pe3 63, (p e (0,7t/2) 

are entangled with negativity N{^^) = sin(p cos(p, but Pas = 
63 (g) 63 is separable. 

When the initial state is the isotropic W, then the nonvanishing 
matrix elements of Was are 

(Was)55 = ^(l+2p), (Was)59 + | 



and 



(Was)66 = ^(l-/'), (Was)88 
(Was)99 = ^ 



(V.4) 



(V.5) 



VI. SYSTEM OF ATOMS OF TYPE (I). NUMERICAL 
RESULTS 



When we consider two atoms of type (I), the analysis of 
evolution of the system is much more involved. For distant 
atoms the master equation reads 



(VLl) 



with 



iYi(2o^iPof3-crfiP-pafi) 

\ y2(20^2p(^3 - <42P - P<42) 

^yi2(2a^iP023 - O21P - PO21) 
^ yi2(2a^2P(^3 - <^2P - P<^2) 
^yi(2ofipof3-ofip-pofi) 



(VI.2) 



i y2(2of2pof3 - of2p - pof2 



yi2(2ofipof3-ofiP 



pofl) 



2'yi2(M2P<3-<2P 



-Po'n) 



where of/, of/ for k,l = 1,2,3 are defined as in Sect. IV. No- 
tice that in the present section, the vectors ei , 62 and 63 rep- 
resent states |2) and |3) respectively and cross damping 
rate y^ is given by ( III.3> . 

From (IVI.U we obtain a system of differential equations 
for matrix elements of density matrix p which we solve nu- 
merically. The results can be used to analyse time evolution 
of negativity of given initial state. Since in a single three-level 
atom the effect of quantum interference can be observed in 
both systems (1) and (11) |9|, we expect that also in this model 
the presence of interference will slow down the process of dis- 
entanglement. 




So for p > i, asymptotic states are still entangled. 



FIG. 6: Time-dependence of negativity for pure initial state <III.3> 
with (p = 7l/6, 9 = n/8, for Y2/Y1 = 0-9 and different values of P/ 
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FIG. 7: Time-dependence of negativity for isotropic initial state 
<III.6> witii p = 3/4, for Y2/Y1 = 0-9 and different values of P/ 



Numerical results entirely confirm these expectations. In 
particular, we show that when P/ < 1, all initial states evolve 
to the asymptotic state |e3)(e3|, so disentangle asymptoti- 
cally, but as in the previous case, the rate of disentangle- 
ment depends on quantum interference: for larger P/ the pro- 
cess of disentanglement is slower FIGS. 6 and 7 represent 
the results of calculations for pure initial state (IIII.3> with 
cp = Jt/6, 9 7t/8 and isotropic state (IIII.6> with p~3/4. 

In the case of maximal interference (p/ = 1), we also ob- 
serve that the system has nontrivial asymptotic states with 
non-zero negativity. For example, FIG. 8 represents the re- 
sults of numerical calculations of evolution of negativity for 
maximally entangled initial state ( IIII.5> . 
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FIG. 8: Time-dependence of negativity for initial state <III.5> in the 
case of maximal interference and for different values of K = Y2/Y1 



VII. CONCLUSIONS 

We have studied dynamical aspects of entanglement 
in a system ot two independent three-level atoms in V- 
configuration, interacting with the vacuum. Spontaneous 
emission causes decoherence and degradation of initial entan- 
glement, but quantum interference between principal transi- 
tion in each atom can slow down the process of disentangle- 
ment. We have shown this effect using two models of dynam- 
ics, which at the level of single three-level atom predict the 
phenomenon of interference. Our results indicate that these 
two kinds of quantum evolutions similarly describe dynami- 
cal behaviour of entanglement in the system. 
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